Preprint typeset in JHEP style - PAPER VERSION HEP-TH/0309 1 10 

Extended QCD 2 from dimensional projection of QCD^ 



Jorge Alfaro 1 , Alexander A. Andrianov 2 , Pedro Labrana 1 

1 Facultad de Fisica, Pontificia Universidad Catolica de Chile, Vicuna Mackena 4860 
Macul, Santiago de Chile, Casilla 306, Chile 

2 V. A. Fock Department of Theoretical Physics, St. Petersburg State University, 1, ul. 
Ulianovskaya, St. Petersburg 198504, Russia and Istituto Nazionale di Fisica Nucleare, 
Sezione di Bologna via Irnerio 46, Bologna 40126, Italy 

E-Mail: jalfaro@puc.cl, andrianov@bo.infn.it and plabrana@puc . cl 



Abstract: We study an extended QCD model in (1 + 1) dimensions obtained from QCD in 
AD by compactifying two spatial dimensions and projecting onto the zero-mode subspace. 
We work out this model in the large N c limit and using light cone gauge but keeping the 
equal-time quantization. This system is found to induce a dynamical mass for transverse 
gluons - adjoint scalars in QCD2, and to undergo a chiral symmetry breaking with the 
full quark propagators yielding non-tachyonic, dynamical quark masses, even in the chiral 
limit. We study quark-antiquark bound states which can be classified in this model by their 
properties under Lorentz transformations inherited from AD. The scalar and pseudoscalar 
sectors of the theory are examined and in the chiral limit a massless ground state for 
pseudoscalars is revealed with a wave function generalizing the so called 't Hooft pion 
solution. 
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1. Introduction 

It is fairly well accepted that Quantum Chromodynamics (QCD) is a correct theory of 
strong interactions giving a good description of the high energy scattering of hadrons. 
However our ability to develop its predictions at low energies, the non-perturbative regime 
of QCD, is not as good. The latter problem can be handled more easily if we consider 
a lower dimensional space. In this way, QCD in two dimensions, QCD2, first studied by 
't Hooft in the large N c limit [1,2], has proved to be a very interesting laboratory for 
checking non-perturbative aspects of QCD. 

This model has been a subject of extensive investigation, covering the study of quark 
confinement, hadron form factors, deep-inelastic [3,4] and hadron-hadron scattering [5], 
the role of instantons and solitons [6,7], multi-flavor states [8], quark condensation [9] and 
heavy quark physics [10]. The bosonization of QCD2 was used as a tool for the description 
of meson and baryon spectra [11-13] and of chiral symmetry breaking [12,14]. 

The remarkable feature of this model is that it exhibits quark confinement with an 
approximately linear Regge trajectory for the masses of qq bound states [2]. As well QCD2 
in the large N c limit has an interpretation in terms of string variables [15] sharing thereby 
important properties that we expect to find for real QCD in four dimensions. 

However QCD2 is too simple to describe certain realistic aspect of the higher dimen- 
sional theory. For instance, there are no physical gluons in two dimensions and both angular 



- 1 - 



momentum and spin are absent in a 2D theory. Moreover in QCD2 there is a problem 
of tachyonic properties of 't Hooft quarks [2] and a realistic chiral symmetry breaking is 
absent which manifests itself in the decoupling of massless pions [16]. Thus one needs 
more information from extra dimensions in order to have a more realistic picture of QCD 
by using a 2D model. For instance, the inclusion of boson matter in the adjoint gauge 
group representation [17, 18] may provide an information of transverse degrees of freedom, 
characteristic of a gauge theory in higher dimensions, and give a more adequate picture of 
strong interaction. 

In this paper we study a QCD reduced model in (1 + 1) dimensions which can be 
formally obtained from QCD in AD by means of a classical dimensional reduction from 
AD to 2D and neglecting heavy K-K (Kaluza-Klein) states. Thus only zero- modes in the 
harmonic expansion in compactified dimensions are retained. As a consequence, we obtain 
a two dimensional model with some resemblances of the real theory in higher dimension, 
that is, in a natural way adding boson matter in the adjoint representation to QCD2. 
The latter fields being scalars in 2D reproduce transverse gluon effects. Furthermore this 
model has a richer spinor structure than just QCD2 giving a better resolution of scalar 
and vector states which can be classified by their properties inherited from AD Lorentz 
transformations. The model is analyzed in the light cone gauge and using large iV c limit. 
The contributions of the extra dimensions are controlled by the radiatively induced masses 
of the scalar gluons as they carry a piece of information of transverse degrees of freedom. We 
consider their masses as large parameters in our approximations yet being much less than 
the first massive K-K excitation. However this model is treated on its own ground without 
any further comparison with possible contributions of heavier K-K boson states. As well 
we do not pay any attention to the mass generation via Hosotani mechanism [19] as it has 
no significance in the planar limit of two-dimensional QCD models under consideration, 
just being saturated by fermion loops sub-leading in the 1/N C expansion. 

This model might give more insights into the chiral symmetry breaking regime of 
QCD 4. Namely, we are going to show that the inclusion of solely lightest K-K boson 
modes catalyze the generation of quark dynamical mass and allows us to overcome the 
problem of tachyonic quarks present in QCD2- 

The paper is written as follows, in Sect. 2 we make a formal design of the 2D reduced 
model from QCD in AD by using Kaluza-Klein compactification and further projection 
onto zero-mode subspace. First we work with the gluon part of QCD4 which is mapped 
into a 2D SU(N) gauge field and two scalar fields in the adjoint representation of the 
reduced gauge model. In principle, these scalars are massless but we recover an infrared 
generation of masses for these scalar gluons. Next we examine the fermion part of QCD4 
and after compactifying and projecting we obtain two spinors minimally coupled to the 2D 
gauge field and interacting with the scalar gluons. The full propagator for these fermions 
is calculated in the large N c limit by using the Schwinger-Dyson equation and a non- 
perturbative solution for the self-energy is found that gives a real net quark mass, even in 
the chiral limit [20]. In Sect. 3 we study bound states of quark and antiquark by using 
the Bethe-Salpeter equation in the large N c and the ladder approximations, qq states of 
the model are identified in terms of their custodial symmetries from AD Lorentz group. 
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We focus our attention on the scalar and pseudoscalar sectors of the theory and in Sect. 
4 the low mass spectrum of the scalar and pseudoscalar sector is examined. In the chiral 
limit, a massless ground state is revealed which is interpreted as the pion of the model, 
a projection of the real Goldstone boson in 4D to 2D treated in the large- A^ c limit. We 
summarize our points and give conclusions in Sect. 5. 

2. Compactification of QCD A to the Reduced Model 
2.1 The Gluon Part 

We start with the QCD4 action for gluons in (3 + 1) dimensions: 

Sqcd = - J d 4 x^tr(Gl u ), (2.1) 

where: 

Gjxv = d^G, - d u G^ -i[G„, G„] , (2.2) 
with G^ = G a ^\ a . Generators A a and structure constants f a bc satisfy: 

tr(\ a \ b ) = ^5 ab , [A , A 6 ] = ifabcK • (2.3) 

Now we proceed to make a dimensional reduction of QCD, at the classical level, from 
4D to 2D. For this we consider the coordinates ^2,3 being compactified in a 2- Torus, 
respectively the fields being periodic on the intervals (0 < £2,3 < L = 2ttR). Next we 
assume L to be small enough in order to get an effective model in 1 + 1 dimensions. 

Following this scheme we expand the fields in Fourier modes on the compactified di- 
mensions: 

00 {2m ] 

G^(x ,x 1 ,x 2 ,x 3 ) = ^2 exp < — [nx 2 + mx 3 ] \ G"' m (x , xi). (2.4) 

n,m=—oo ' 

By integrating heavy K-K modes Eq. (2.4) with n,m / in the QCD generating 
functional we obtain the low-energy effective action for zero-modes which can be efficiently 
prepared along the Wilson construction with the help of a cutoff A. Accordingly, let us 
adopt the low-energy limit |po,i| <C = ^Vn 2 + m 2 ; n,m/0. Hence if we choose 

a 2D cutoff A ~ i = ^ larger than the typical momenta po : i, we are led to analyze the 
zero Fourier modes only. Thus all relevant fields on the uncompactified 1 + 1 space-time 
consist of zero modes. 

To make contact with 't Hooft solution, we introduce light-cone coordinates in the 
uncompactified space-time x± = "^(^0 ± xi), then we obtain A± = -^(Go ± G\) and two 
scalar fields defined as (f)i )2 = G 2 ^. 

By keeping only the zero K-K modes , we get the following effective action: 

r t 2 

S= d 2 x— tr[(d-A + -d+A„ -i[A_,A + ]) 2 + 2(d-4> 1 - i[A_,</, 1 ]){d + ^ 1 -i[A+,<f> 1 ])+ 
J 9 
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2(d-<h-i[A-,<h])(d+h-i[A + ,<h]) + [<h,<h?] • (2.5) 

Following the t'Hooft scheme we choose the light-cone gauge = 0. Then the action 
is transformed to: 

r j2 

S= / d 2 x— tr[(d^A + ) 2 -2iA + ([4> 1 ,d^ 1 ] + [4>2,d^ 2 ]) 
J 9 

+2(d^id+cpi + d-fad+fa) + [0i,0 2 ] 2 ] • (2.6) 

This action is quadratic in A+. To analyze the mass spectra for fields 01,2, we integrate 
over A + : 



f L 2 

S 4> = J <? x -=ptr{2d-(t> k d + <j) k + fc] 



I 



d 2 y {[fa, d-<h] + [02, d-H) (x\P.V.^\y) (fad-fa)}. (2.7) 



We expect the infrared mass generation for the two-dimensional scalar gluons 0i ;2 [21]. 
Let us calculate it to the leading order of perturbation theory. Applying the Background 
Field Method we split = j^(4> c + ?7), with C being a constant background. Then the 
one-loop contribution is given by: 

5( 2 ) = J d 2 xtv[2d- Vk d +Vk + ^([hcVi] 2 + [<t>icM 2 ) (2-8) 

g 2 f 1 

~j2 j d 2 yd^([(f) kc ,r] k (x)})(x\-^\y}d-([(l) kc ,ri k {y)]) 



f 1 

+ 2 ^2 (ivumi mc 4>2c] + [(f>ic,m][m,<l>2c])f- 



The masses can be calculated from the quadratic part of the effective potential, after 
integration by parts: 



5( 2 ) = J d 2 xtr{2d-r n d + r n + 2d^ 2 d+m 



7,2 



+j2{[<f>ic,m} 2 + [02c, vi} 2 + [0ic V2] 2 + [02c, m} 2 )}, (2.9) 

They receive equal pieces from the two-dimensional gluon exchange and the tadpole. 

To estimate the masses of scalar gluons we use the Schwinger-Dyson equations as 
self-consistency conditions: 

r2 2N c g 2 f A d 2 p 1 N c g 2 A 2 T _A 2 + M 2 



thus M 2 brings an infrared cutoff as expected. We notice that the gluon mass remains 
finite in the large- N c limit if the QCD coupling constant decreases as 1/N C in line with 
the perturbative law of AD QCD. Respectively such a mass would dominate over the 
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mass induced by the Hosotani mechanism as the latter one is sub- leading in 1/N C if being 
provided by fermions in the fundamental representation [19]. 

We adopt the approximation M^;A~l/i?to protect the low-energy sector of the 
model and consider the momenta |po,i| ~ M. Thereby we retain only leading terms in 
the expansion in and and also neglect the effects of the heavy K-K modes in the 
low-energy Wilson action. 

Now we proceed to the evaluation of the model (2.5) as a generalization of the 't Hooft 
QCD2 and take a bounded value of M. Then the low-energy approximation corresponds 
to the small coupling constant: 



M 2 ' 

c iU & TP 

in accordance to Eq. (2.10). Thus the relevant coupling constant for the model (2.5) is: 

, N c f- N C ?A* M 2 

This dimensional constant governs the meson mass spectrum of the 't Hooft QCD2 [1] and 
can be associated to a particular meson mass, i.e. be fixed. Then the gluon mass is slowly 
growing with the increasing cutoff A and and the "heavy-scalar" expansion parameter is: 

A = 9 , r0 = — 1 -tt. (2.13) 

We observe that the limit Mc A supports consistently both the fast decoupling of the 
heavy K-K modes and moderate decoupling of scalar gluons, the latter giving an effective 
four-fermion interaction different from [23]. Nevertheless we prefer to keep first the gluon 
mass finite and further on derive the expansion in parameter A gradually from the Bethe- 
Salpeter equation because for small A one recovers non-analytical infra-red effects for the 
ground state wave function. 

2.2 Introducing fermions 

Now let us make the K-K reduction of the quark part of QCD4 at the classical level. The 
original Lagrangian reads, 



/ 



d 4 x$> (i'f - m) * , (2.14) 

where the covariant derivative is defined conventionally: = — i . The following 
representation of 7 matrices is employed: 



<ri w 2 / ia 3 \ 

7o = n , 7i = n . , 72 = . „ ,73 





suitable for the purposes of K-K reduction. Then the bispinor is decomposed into two- 



dimensional ones, ^ = [ a ] with the help of the projectors, P a j = ^(1 ± 7275). 
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Figure 1: The quark propagator in terms of the one-particle irreducible self-energy where it is 
defined £ = i\- 

Now we compactify on the torus and expand the fermion fields in Fourier modes of the 
compactified dimensions £2,3: 

^(x ,x 1 ,x 2 ,x 3 ) = Y] e^pl-^-(nx 2 + mx 3 )\^ m (x ,x 1 ). (2.15) 



n,m=— 00 



As well as for gluon fields, we retain only the zero modes, obtaining two spinors ip a 
and ipb of two components each for the model in (1 + 1) dimensions. 

By keeping the zero K-K modes only, one arrives to the following effective Lagrangian 
for the quark part of the model in (1 + 1) dimensions, after a suitable rescaling of fermion 
fields: 



£ - ir„ ( -'Kdn-i-^j" A,, } ('„ - mi;, r„ + (a - b) + 

9 



4>1 {4>a (io-3) 4>b + A (^3) Ipa) + 7=== 02 (Va (^3) ^ a - (iff 3 ) , (2.16) 



where we have used the definition for the 7 matrices in two dimension: 7^ = (7°,7 1 ), 
7 = 01, 7 1 = ici2, inherited from (2.15). 

The 7 matrices in the light cone basis are: 7 =t = 4= (7 ± 7 1 ) and satisfy {7 + , 7~} = 2 

and ( 7 +) 2 = ( 7 -) 2 = 0. 

We consider the dressed quark propagator, which can be expressed in term of one- 
particle irreducible self-energy \ as follows (see fig. 1) : 

Sip) = ? " -^TT • (2-17) 

[7- P+ + 7+ P- ~ m - tx(Pl\ 

The equation for \ ls given graphically in fig. 2, in the large N c limit and in the one- 
boson exchange approximation. To solve the equation we decompose the matrix \'- 

X = X+1- + X-7+ + XoH, (2.18) 

where x+>X-> anci Xo ar e functions of p. The equations for the components (2.18) of \ 
read: 

= 2 2 f [q-ixL p_ (2 19) 

X+ J (2vr)2 2 [q - i x } + [q - i x }_ - (m + i X o) 2 + ie (p - ?) 2 



+ 2 5 2 



d 2 <? 1 [?-*x]- 



(2tt) 2 [2(p - g)+(p - g)_ - M 2 + ie] 2 [ g - i X }+ [q - »x]_ - (m + ixo) 2 + ie ' 
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X- = 2g 2 



d 2 q 1 [q-ix]- 



Xo = 2g 2 



(2vr) 2 [2(p - q)+(jp - g)_ - M 2 + ie] 2 [g - i X ]+ [g - ix]_ - (m + i X of + ie ' 

(2.20) 



-1 [m + ixo] 



(2vr)2 [2(p - q) + (p - g)_ — M 2 + ie] 2 [g - i X } + [q - i X L - (m + i X o) 2 + ie ' 

(2.21) 



2 



We intent to solve these equation for \p\ <C M and at first sight one finds X+ = f~i 
X- = 0, and xo = as solutions, that is, a 't Hooft solution for x+ an d vanishing new 
terms. However it happens that this solution entails imaginary masses for some meson 
states if one takes into account the scalar field exchange, see appendix A. 

On the other hand there exists another possibility, namely, if one looks at Eq.(2.21) 
one could grasp that when the current quark mass m = 0, still a non-zero real solution 
for So = ixo exists, even in the case when M 2 — > oo. To explore this possibility, we start 
with X- = and a constant xo ^ 0, then in the zero approximation in expansion one 
arrives at the 't Hooft solution for x+- 

. 2 f dq- P i g 2 
X+ = -mg / -—2 sgn{q-)- -y = . (2.22) 



(2vr) 2 * KH -'{p-q) 2 _ 2vr p 

To obtain xo m the chiral limit we treat Eq.(2.21) in the low-momentum limit (M 2 3> 
p 2 ), then the following equation determines xo'- 

2 f d 2 q -i 1 _ (f_ 1 (M 2 \ 

9 J (2T) 2 [g 2 — M 2 + ie] [g 2 - m 2 + ie] ~~ 2tt M 2 - m 2 ^^j ' ( } 



where we have defined m 2 d = So — . If one assumes M 2 » m 2 , one can solve the equation 
analytically to find Sq: 



■ 9 2 9 , 2tt 2^ 9 2 M A g 2 
S 2 , ~ — + M 2 exp T M 2 ) ~ — + — 

7T C? 2 7T A 2 



(2.24) 

7T 



Let us notice that in virtue of Eq.(2.12) in this limit So does not tend to zero but 

2 

becomes a constant ^ which fully compensates the tachyonic quark mass of the t'Hooft 
solution. We see also that the solution Eq.(2.24) fully justifies the assumption M 2 » m 2 ,. 
Thus in our approximation m 2 , = in the chiral limit. 

The quark propagator on this solution takes the following structure: 



„, , _ i (7- [P ~ ix] + + 7+ \P ~ ix]_ + m + ixo) J- 
" \P) ~ ~~ 7^7 ^ — 7 7 1 — : v?~, — : — — 1 



+ 7+ p- + S 



2(p- «x)+(P - *X)- - ("i + iXo) 2 + ie p 2 - m 2 , + ie 

(2.25) 
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Figure 2: The equation for E. 



Notice that the dynamical mass in the propagator is real and we don't have tachyonic 
quarks as in [2]. The existence of tachyonic quarks had been interpreted as a signal of 
confinement, but this is incorrect as it was first understood in [3,4,22]. 

3. Mesons 

3.1 The Bethe-Salpeter equation 

Now we want to study qq bound states. In our reduction we have four possible combi- 
nations of quark bilinears to describe these states with valence quarks: (ipaipa), (ipa^b), 
{ipb ^a); {ipb i>b)- Each of these bilinears can be, of course, supplemented with a number 
of gluons and adjoint scalars. Moreover, unlike the t'Hooft model, the interaction in our 
dimensionally reduced model mixes all sectors of the Fock space with quark, antiquark 
and any number of scalars, even in the planar limit, as it has been realized before in the 
light-cone quantization approach [18]. However in the low-energy region for M 2 > 
one can successfully analyze the deviation from the t'Hooft bound state wave function just 
retaining the ladder exchange both by gluon and by scalars for calculation of the valence 
quark wave function and neglecting vertex and self-energy corrections as well as the mixing 
with hybrid states which are (superficially) estimated to be of higher order in the 1/M 2 
expansion. In this approach the qq wave functions satisfy the Bethe-Salpeter equation, for 
example ip a ^ a satisfies an equation given graphically in fig.3. The Bethe-Salpeter equa- 
tions for these four states mix them (see fig.3) but if we define the following combinations 
the equations decouple: 

HP, q) = \ {^ab - ^ba) , A(p, q) = ^ + V bb ) , (3.1) 

9(p, q) = - (^ab + *6a) , 0(p, q) = 2 " ®bb) ■ (3.2) 

These combinations have a direct interpretation in term of Dirac bilinears of the theory 
in (3 + 1) dimensions. In particular, the scalar and pseudoscalar densities, which are within 
our scope in this paper, can be represented by, 



S = ^ = 2Tr(A), (3.3) 
P = ^ 7 5 ^ = -2iTr(Aa 3 ) , (3.4) 



were ^ is a four dimensional bispinor. 
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Figure 3: Bethe-Salpeter equation for the state ^ oa . Notice that in this equation the state 
appears. 



Respectively the manifest form of the Bethe-Salpeter equations on wave functions (3.1) 



is: 



A(p,q) = iS(p-q) (igj^ 



d 2 k 



(k - qf_ 



Hp, k) 



(igj-) S(-q) 



2iS(p- q) (ga 3 



d 2 k 



[2(p-q) + (p-q)_ -M 2 + ie] 



a(p, k) 



(g<T3)S(-q). (3.5) 



A(p,q) =iS(p-q) (igj-) 



d 2 k 
(2vr) 2 



Q-qf_ AM 



(igi-)S(-q) 



+ 2iS{p-q) (ga 3 ) J 



d 2 k 
J2k) 2 



[2{p - q) + {p - q)- - M 2 + ie] 



A(p,k) 



(9*3)S(-q). (3.6) 



Further on we are going to focus our analysis on these equations to describe scalar and 
pseudoscalar meson states. 

3.2 Integral Equation for A and Pseudoscalar Mesons 

Let us explore solutions of the matrix integral equation (3.5). To do that we decompose 
the matrix A as follows: 



HP, Q) = A+ 7+ + A_ 7_ + A + _ 7+ 7_ + A_+ 7_ 7+ , 
then in components Eq.(3.5) yields four integral equations: 



A+(p, q) 



-^g 2 {p-q)-q- 



D{p-q}D[q] J (2vr) 2 



d 2 k 



(k ~ q)i 



A+(p, k) 



(3.7) 



(3.8) 



+ 



d 2 k 



2(p - g)_<Z_A_ - 2g- S A + _ + 2 S (p - 9 )_ A_. 



S n 2 A, 



D[p-q]D[q]J (2vr) 2 



A-(p,q) 



[(k — q) 2 — M 2 + ie] 



2ig 2 S 2 



d 2 k 



D[p-q]D[q\J (2^) 2 



(k-q) 2 _ 



A+(P, k) 



(3. 
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+ 



2ig 2 



D\p-q]D[q}J (2tt) 2 [ 



d 2 k 



2[p - q] [q]A + - 2[q] S A_+ + 2 S [p - g ]A^ 



S n 2 A_ 



[(k - q) 2 — M 2 + ie] 



A+_(p, g) 



2i g 2 S (p - g)_ r d 2 k 
D\p-q}D[q] J (2^ 



A+(p, fc) 



(3.10) 



+ 



2^ 2 



d 2 fc 



D[p-9]£»M7 (2vr) 2 



2(p - g )_ [-g]A_ + + jg] Sp A + - Sp (p - g)_ A_ + S 2 A + . 
[(fc - g) 2 - M 2 + ie] 



2ig 2 Sp g_ /• d 2 fc 
-+^ q )~ D\p-q]D[q]J (2vr) 2 |_(fc-g)? 



A+(p, k) 



(3.11) 



2^ 2 



d 2 fc 



-2[p - g]g-A + - + g_ S A_ - S jg - q]A + + S 2 A_ + 
[(k - q) 2 — M 2 + ie] 



D[p-q]D[q]J (2vr) 2 [ 

where we have defined D[q] = q 2 — m 2 d + ie, and [q] = q + + |^ (j^z) • 

In order to have control over the expansion in let us introduce the following vari- 



ables: 



A_ (p, g) = -2 (P ^ g A_ (p, g) + A+ (p, g) , A+„ (p, g) = -2^2 A+_ (p, g) + A+ (p, q) , 
A_ + (p, q) = 2 {P ~f- A_ + (p, g) + A+(p, g) . (3.12) 

One can see that in the large M 2 limit the equations for these new variables contain 
only terms, superficially proportional to Meantime the equation for A + contains 

't Hooft operator (which is of zero order in as well as the terms of order in which 
both the new variables and also A + appear. Since the corrections to the 't Hooft equation 
for A + are examined here at the leading order in expansion, we retain in this equation 
only the terms including A + : 



a / \ o- 2 (p- q)-q- 

A +^ =2 ' S D\p-q]D[q] 



r 2 



d 2 k 1 A + (k) S 2 , f d 2 k A+(k) 



(2vr) 2 (p - fc)_ \{k - q) 2 — M 2 + ie] (p - g)_ g_ J (2vr) 2 [(k - g) 2 - M 2 + ie] 



S 2 f d 2 k 1 A+(fc) T% [ d 2 k 1 A+(k) 



(P - ?)- 7 (2vr) 2 fc_ [(fc - g) 2 - M 2 + ie] g_ 7 (2vr) 2 (p - jfc)_ [(it - g) 2 — M 2 + ie 

In the region (k — q) 2 <C M 2 one can integrate over g + on both sides of equation (3.13), 
obtaining the following integral equation for the eigenvalues p 2 : 
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- A L dy ¥^ + A [(i- L dym + A ^l dy [(i • 

(3.14) 

Where we have conventionally introduced 4>{q~) = / dq + A + and dimensionless vari- 
ables g_ = = yp- as well as the expansion parameter A, previously defined in 
Eq.(2.13). The integral equation (3.14) has a symmetric kernel that guarantees that the 
eigenvalues p 2 are real. It is also symmetric under the transformation x — > (1 — x) found 
by 't Hooft in [2]. 

To the first order in A the straightforward calculation from Eq.(3.13) gives = 0, 
but such a kernel makes this equation too singular to have any reasonable solutions. This 
singularity must be softened by re-summation of a series of similar singularities of higher 
order in A. Just in order to simulate this re-summation we have introduced (3 which 
effectively regularize the integral equation. The coefficient (5 governs the behavior of the 
solution of Eq.(3.14) near the end points. We fix it by imposing self consistency of solutions 
near the boundaries. Although the parameter (3 is small and vanishing with ^4^0, its 
presence guarantees the positivity of the spectrum (see appendix B). 

4. Integral Equation and Approximate Spectrum 

To explore solutions to Eq.(3.14) we examine small A, which allows us to put to 2 , = 0. 
Working in this limit we select out consistently the wave functions that behave as x 13 or 
(1 — x)P on the boundaries, then (3 is fixed to in order to cancel the principal divergences, 
x^^ 1 and (1 — x)P~ 1 , that appear when x is near the end points. 

Evidently Eq.(3.14) does not mix even and odd functions with respect to the symmetry 
x -4=>- 1 — x. On the other hand the ground state should be an even function. When 
inspecting the wave function end-point asymptotics from the integral equation (3.14) one 
derives the following even function as a ground state solution for A — > limit: 

<fo(z) = (4x[l - s])t - - (4x[l - x])? . (4.1) 

7T 

This is basically a non-perturbative result that differs from 't Hooft solution in the 
A — > limit, giving p 2 = 0, as we would expect from spontaneous chiral symmetry breaking 
in 4D. Notice that the second term in Eq.(4.1) is a leading contribution to the ground state 
solution for 't Hooft equation with to 2 = ^ [2], and the first term compensates powers of 
A when we insert 4>o{x) in Eq.(3.14) and take the A — > limit. The entire result is a shift 

2 2 

of the "unperturbed" ground state p 2 = 7.25^7 G us t the 't Hooft solution with to 2 = 
to the non-perturbative solution with p 2 = in the A — > limit. For the other massive 
states we are unable to find analytic solutions, as happens with 't Hooft equation, but we 
could estimate them working with the Hamiltonian matrix elements p 2 ((p,(p) = (cp,H(j)) 



to; 



p - 



x 



m„ 
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and using the regular perturbation theory, starting from 't Hooft solutions (A = 0) [2]. 
Formally the set of solutions to Eq.(3.14) is composed of even and odd functions with 
respect to the reflection x -4=>- 1 — x, but only the even functions have physical meaning 
as describing pseudoscalar states. It occurs because the 4D pseudoscalar density could be 
written as a combination of 2D pseudoscalar densities for which it is possible to use the 
arguments given in [3]. Thus for the pseudoscalar sector only the set of even solutions of 
Eq.(3.14) is selected out. Similarly for scalar states the equation differs from Eq.(3.14) by 
signs of perturbation terms, i.e. A — ► —A and the appropriate set of solutions consists 
of odd functions, then only perturbative massive solutions are possible for these states. 
As we have significantly departed from QCD4 neglecting heavy K-K states and taking the 
limit A — ► we don't make any numerical fits to some observable mesons with quantum 
numbers of scalar and pseudoscalar radial excitations just postponing this work till a more 
systematic estimations of neglected terms. 

5. Conclusions 

We summarize our points. Quantum Chromodynamics at low energies has been decom- 
posed by means of dimensional reduction from (3 + 1) to (1 + 1) dimensions and a low 
energy effective model in 2D has been derived. It is presumably more realistic than just 
taking QCD in 2D because the model includes, in a nontrivial way, physical (transverse) 
gluon degrees of freedom. We treat this model in the equal-time quantization approach 
using Schwinger-Dyson equations instead of Hamiltonian quantization, and we argued that 
the adjoint scalar fields dynamically gain masses. We did an explicit analysis of the meson 
bound states by using Bethe-Salpeter equation in the limit of large scalar masses. 

We found that the model has two regimes classified by the solutions of quark self- 
energy. In one case, the perturbative one, we have tachyonic quarks in the chiral limit 
as in [2] and our model could be interpreted as a perturbation from the result of chiral 
QCD in 2D. But we are not allowed to consider that possibility because the spectrum for 
the lowest qq bound states becomes imaginary once we consider small but nonzero A, see 
appendix A. This means that the 't Hooft chiral solution is not a stable background to 
start a perturbation to QCD2 that includes transverse degrees of freedom. Also in this 
regime the equations that describe the mass spectrum for scalars and pseudoscalars are 
equal, certainly we have oversimplified something. On the other hand the non-perturbative 
solution to quark self-energy supports non-tachyonic quarks with masses going to zero, in 
the chiral limit. 

We focused our attention on bound states of quark-antiquark in the valence quark 
approximation which is right for low energies. It was shown that qq states constructed in 
our reduced model have a direct interpretation in terms of their properties under Lorentz 
transformations inherited from (3+1) dimensions. Respectively, we have different equations 
that govern the scalar and pseudoscalar mass spectrum and in both cases the spectrum is 
positive definite with a massive ground state for scalars and a massless ground state for 
pseudoscalars. 
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The pseudoscalar and scalar sectors of the theory have been analyzed and a massless 
solution for the pseudoscalar ground state has been found. We interpreted this solution as a 
"pion" of our model although it is a quantum mechanical state only existing in the large- N c 
limit. The appearance of massless boson in a two dimensional theory in a particular limit 
is not a new feature and has been studied in [9, 13,25] and more recently justified in [26]. In 
particular in the planar QCD2 the 't Hooft pion solution is a massless boson bound state, 
related with a non-anomalous U(l) symmetry. We expect that our pion solution Eq.(4.1) 
has more physical properties from the AD theory, as compared to the massless solution of 
QCD 2D. For the latter one, it has been proved that the 't Hooft pion state is decoupled 
from the rest of the states [16]. On the contrary, the argument of non-anomalous U(l) 
symmetry that produced the decoupling of the massless scalar field in QCD2 does not hold 
in our model. 

Indeed in QCD2 the U(l) vector current = ■ip a j fJ "ip a + ■ip b 'j fl ip b is conserved then 
= 'Q v (j)(x). Due to 7^75 = t l ivl v we can write the U(l) axial current as follows 
Jjj, = <9 M </>(x), then if the axial current is also conserved we have d^d^cfi = 0, a massless 
decoupled scalar field. This scalar field is identified with the 4>o = 1 't Hooft state. 

In our model the situation is different because the U(l) axial symmetry is explicitly 
broken by the interaction term in Eq.(2.16), but our model has another symmetry due 
to the breakdown of the AD Lorentz symmetry ,50(1,3) to 50(1, 1) x 0(2), where 0(2) 
now becomes an internal symmetry. The vector and axial currents associated with this 
symmetry are conserved and the axial symmetry corresponds precisely to the projection 
of the U(l) axial symmetry in AD to our model in 2D. Our massless pseudoscalar meson 
is related with these currents, then we cannot expect the same behavior for our massless 
pseudoscalar as in QCD2 because we are talking about different symmetries. For instance, 
the vector and axial 0(2) currents do not satisfy the same relation as and in QCD2. 
On the other hand the possibility to have a coupled massless scalar in the large--/V c does 
not contradict the Coleman theorem as it has been first pointed out in [25], and particular 
in our case because it is not a quantum Goldstone field but a quantum-mechanical state. 
We stress that this composite massless state can be related to the real pion, to some extent, 
at the tree level as the chiral 2D perturbation theory being renormalizable is drastically 
different from the QCD AD Ch.P.T. [27]. 
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A. Quark-Self Energy, First Solution 



Let us consider the t'Hooft solution for the quark self energy \+ = 2^|~> Xo = an d 
X- = and explore what this solution implies for the bound state spectrum of quarks 
and antiquarks. We proceed in the same way as before obtaining a similar Bethe-Salpeter 
(BS) equation for A(p,q) Eq.(3.5), but now the dynamic mass is the one found in [2] 



m 



— ) and Eq is zero. In the chiral limit 



m 



0), we write the BS equation in 



terms of the decomposition given in Eq.(3.7) and obtain two decoupled sets of equations: 
one for A + (p,q) and A_(p, q) and the other one for A + _(p, q) and A_ + (p, q). To calculate 
corrections to the 't Hooft equation for A + we retain just the first set of equations: 



A+(p, q) 



-4ig 2 (p- 



D[p - q] D[q] 



(A.l) 



111 
J (2vr 



(2vr) 2 (k - q)l 



A. 



.(p,k)-J 



d 2 k 



A-(p,k) 



(2vr) 2 [(k - q) 2 — M 2 + ie] 



A-(p,q) 



Aig 2 [p — q][q] f d 2 k 



A + (p,k) 



D\p - q] D[q] J (2vr) 2 [(k - q) 2 — M 2 + ie] ' 
where we can write [q] = \{q 2 — m 2 ,) {j^j- Then the equation for A_(p, q) becomes: 



(A.2) 



A_(p,fe) =ig 2 



ip-k). 



2U 



A + (p,k') 



(2vr) 2 [(k - k') 2 — M 2 + ie] ' 



(A.3) 



now we replace this equation into Eq.(A.l) and integrate over k in the second term retaining 

2 

only the first order. Further on, it is possible to integrate over q + at both sides of the 
equation obtaining the following integral equation for 4>{x): 



P 2 + 



x(l — x) 



dy 



(x - y)< 



<Kv) 



pi p 
+ A / dy- 

Jo {x - y) 



In 



Q--x)y 
(l-y)x 



(A.4) 



where we have defined p 2 = -^iV 2 - O ne might solve this equation by using perturbation 
theory near 't Hooft solution (A = 0). We know from [2] that the ground state of 't Hooft 
equation is (f>o{x) = 1 with p 2 = 0, but the perturbation term in Eq.(A.4) generates a 
negative correction to p 2 . Thereby for non-zero values of A we obtain for the ground state 
an imaginary mass. 
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B. Positivity of the Spectrum 



Here we show that the spectrum of integral equation (3.14) is indeed positive as stated 
above. We know that the lowest state should be an even function. For any function that 

A A 

goes to zero faster than or (1 — x)~^ when x goes to zero or one, respectively, we can 
consider the terms proportional to A as perturbations to 't Hooft term which is positive 
definite. Then the only possible trial functions that could generate negative p 2 are nearly 
constant functions. Let us explore this possibility and consider a function f(x) = [x(l— x)] 7 , 
with 7 small but arbitrary, then we multiply both sides of Eq.(3.14) by f(x) and integrate 
over x obtaining (see appendix of reference [24]): 



p 2 (f(x)\f(x)) = 2^ 



-2 



TM 4 
r[2 7 p 



2Air 



rb + /3] 2 r[ 7 + /?]r[i + 7 ] 



If we assume A and 7 = aA to be small, only the first two terms of the right hand side of 
Eq.(B.l) are relevant. The divergent part of Eq.(B.l) is: 



p 2 {f(x)\f(x)) 



1 

7 



1 - 



2Aj 



(7 + /3) 5 



1 

7 



1 



2a 



l\2 



(B.2) 



Then for any value of a the expression between brackets is evidently positive. 

Notice that the existence of a non-zero (3 guarantees the positivity of the spectrum, 
when 7 goes to zero and A is small but different from zero. 
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